

* counts from counterfactual *
sum pcount if ltd_iss_all_other > 9500 & ltd_iss_all_other <= 10000 
scalar cfbcount = r(sum)
sum pcount if ltd_iss_all_other > 10000 & ltd_iss_all_other <= 14500 
scalar cfacount  = r(sum)

* actual counts  * 
sum count if ltd_iss_all_other > 9500 & ltd_iss_all_other <= 10000 
scalar actbcount = r(sum)
sum count if ltd_iss_all_other > 10000 & ltd_iss_all_other <= 14500 
scalar actacount = r(sum)

* coefficients - Excess Mass
local totmanip_ex = (actbcount - cfbcount) / totalcount 
local inranmanip_ex = (actbcount - cfbcount) / cfacount

* coefficients - Missing Mass
local totmanip = (cfacount - actacount) / totalcount 
local inranmanip = (cfacount - actacount) / cfacount